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TECHNICAL NOTE 3146 


NOTE ON THE AERODYNAMIC HEATING OF AN OSCILLATING SURFACE 


By Simon Ostrach 


SUMMARY 


An analysis of the temperature distributions in a fluid over an 
oscillating surface with heat transfer is made and associated heat- 
transfer parameters are compared with those for the case of conduction 
at a stationary surface with the same initial temperature potential. 
It is found that the heat transfer for the oscillating surface can be 
considerably different from that for conduction alone. The effect of 
the surface oscillations on the thermal state of the fluid is studied 
by means of average static- or total-temperature defects, and it is 
demonstrated that the oscillations could alter the fluid temperature 
appreciably. 


INTRODUCTION 


The increased improvement of present-day propulsion systems and 
the development of new propulsion systems have posed numerous new prob- 
lems in the field of heat transfer. Elucidation of unusually high heat- 
transfer coefficients which sre apparently encountered in unsteady flows 
and means of increasing heat-transfer coefficients under given condi- 
tions are greatly desired. As a preliminary attempt to gain insight 
into such problems, it seems worth while to consider the heat-transfer 
aspects of the classical problem wherein the fluid motion is induced by 
oscillating a conducting surface axially in viscous fluid. In particu- 
lar, the effect of disturbing the equilibrium (steady state) conditions 
after the periodic motion of the fluid has been established will be 
studied. The temperature distributions in the fluid are determined as 
exact closed-form solutions of the energy equation pertinent to the 
problem and, hence, related heat-transfer parameters can be compared 
with those for a stationary surface to demonstrate the effect of the 
surface oscillations. Other exact solutions of the energy equation for 
somewhat analogous problems are presented in references l to 5. 
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ANALYSIS 


The laminar motion of an incompressible viscous fluid induced by 
the axial (longitudinal) oscillations of a bounding plane surface doubly 
infinite in extent is well known and is reported in numerous sources 
(see Schlichting, ref. 4, for example). The equations of motion for 
this case, assuming constent physical renee reduce to the class- 
ical heat-conduction equation 


ou dfu 

E (1) 
ot Sy2 

where u is the velocity component parallel to the surface, t denotes 

the time, v is the kinematic viscosity, and y is the coordinate nor- 

mal to the surface. (All symbols are defined in the appendix, ) The 

associated boundary conditions are - 


u(O,t) = U cos nt (2) 
and. 
ufs, t) = O (3) 


where U denotes the amplitude and n the frequency of the oscilla- 
tions. The velocity distribution, after sufficient time has elapsed for 
the periodic motion to be established, is given by 


uly,t) = U exp (- -/n/2v y) cos (nt - x/29 y) (4) 


The appropriate energy equation is l 
a 
sit Ft (S) (5) 


where T denotes static temperature, a is the thermal diffusivity, and 
is the specific heat at constant pressure. Note that the last term 


on the right of equation (5) is due to aerodynamic heating and describes 
the increase in fluid temperature caused by the oscillations. In order 

for the problem to be physically reasonable, however, after a long per- 

iod of time the system should be in a state of thermal equilibrium (that 
is, in essentially a steady state) and it is thus necessary to cool the 

surface. The thermal equilibrium state can be described by first solv- 

ing equation (5) subject to the boundary conditions 


T(0,t) = Ty (6) 
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and 
T(e,t) = T. (7) 


that is, the oscillating surface is maintained at some uniform tempera- 
ture T,. It is unnecessary here to specify any initial conditions, for 
they will not influence the results after a sufficiently long time. 


After the explicit form of the nonhomogeneous or aerodynamic heat- 
ing term in equation (5) is detexmined by means of equation (4), a par- 
ticular solution of equation (5) was found to be 


ny =~ LPF exp ( Yeah x) fa + ee = VB a) (8) 


4e = Pr 
p 2 


where Pr = v/a is the Prandtl number. The complementary function to 
be added to equation (8} in order that the camplete solution satisfies 
the boundary conditions given by equations (6) and (7) is 


- Pr 
3 2 


T = T, +02 1 + SPL Vale y) og fant - nfa y) - ext (7 


+ (Ta - I) erf (z) (9) 


where erf( ) is the error function and is defined as 


r B 
erf(p) = "Ai expl-rf} dr 
O 


so that 


erf(0) = O 


erf(») = 1 


The complementary function is synthesized from well-known solutions of 
the heat equation; for example, the term containing the product of the 
exponential and trigonometric functions corresponds to the solution of 
the heat-conduction equation for the case where the temperature on a 
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stationary surface is specified to be a harmonic function of time (see 
ref. 5, p. 47), and the two error-function terms combined correspond to 
the solution of the problem of unsteady heat conduction to or from a 
stationary surface due to a temperature difference or potential of 

a = 

a - ge - Tyr Note that the error-function terms describe the tran- 
sient heat transfer which will vanish for large time. The temperature 
distribution in the equilibrium state (that is, after a long time) should 
be independent of time except for periodic terms and is given by the sum 
of equations (8) and (9) with the error-function terms vanishing as 


eof mF) Eo Ay 
(AE Jole AE) ” 


The relation for thermal equilibrium between the surface temperature Te 
and the ambient temperature T, can be found by evaluating equation 
(10) as y> œ and is given by 


T = P ao u“pr (11) 
w eo 4e 
P 
Thus, as was anticipated, it can be seen from equation (11) that for 
thermal equilibrium the surface must be cooled to compensate for the 


aerodynamic heating, and the temperature distribution above a surface 
oscillating for a long time in a fluid at temperature TY is given by 


TaT = pa eno( E y) s = ESE y)ooe( an = 1E v) 
_ =o- E y)eos (ant _ 2 y) (12) 


Now that the fully developed state of the fluid is completely de- 
scribed (by eqs. (4) and (12)), consideration can be given to the prob- 
lem at hand: namely, the effect of the oscillations on the heat trans- 
fer if, at some time (say t = 0) after the periodic motion (steady 
state) is established, the surface temperature is fixed at some temper- 
ature T. different from that given by equation (11). The problem to be 
solved to answer this question is mathematically identical to the one 
for thermal equilibrium except that Tọ, replaces TI, in the boundary 


conditions and that the solution of the present problem must match that 
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given by equation (12) for t = 0; that is, the proper initial condi- 
tion must be satisfied by the temperature distribution in order to den- 
onstrate the transient effects properly. It was previously shown that 
the sum of equations (8) and (9 


pea, + E 1 - em (- fE v) - exe (35) 


+ (To - a Jere MA (13) 


satisfied the differential equation (eq. (5)) and the boundary condi- 
tions (eqs. (6) and (7)), where T, is now replaced by Ty- Evalua- 
tion of equations (12) and (13) at t=O shows that they are identi- 
cal so that the proper initial condition is satisfied and thus equation 
(13) represents the temperature distribution over an oscillating sur- 
face which is maintained at a constant temperature Ti; Because E 

is a surface temperature different from that for thermal equilibrium, 
the solution given by equation (13) is of greatest physical significance 
only as long as the transient terms are important. This of course, is 
the problem of primary interest herein. 


Although the solution given by equation (13) was developed for an 
incompressible viscous fluid with constant property values, it is equally 
valid in the case of a compressible viscous fluid if the boundary-layer 
assumptions are made, if the Prandtl number and the product of p and 
p are taken to be constants, and if y is replaced by 7 where 


"E 
7 = pi p. dy 
O 


This follows because the von Mises transformation (see ref. 6, for 
example) under these assumptions reduces the compressible boundary- 
layer equations to the forms of equations (1) and (5). 


Any pertinent heat-transfer quantities can be obtained from equa- 
tion (33). However, since the time dependence itself is usually not of 
primary practical importance, the temperature gradient at the wall 
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(which is a measure of the heat transfer) will be averaged over a period. 


This average gradient is given by 


ox 
eens tot 


(Z) = 


PeT) to t2. Vi (14) 


K ita 


Ky 
AAT 
AY 
Sep 
= 
YB 
4 
J 
J 
ct 

O 
+b 
aa EY 
i 
I 
S] 

O 


to 


where to indicates the beginning of any cycle. For the equilibrium 
state, an analogous equation is obtained from equation (12): 


oTe fpr Ta 
(Zs) - EE (s) 


This average temperature gradient is due entirely to the aerodynamic 
heating and is also essentially (that is, except for the multiplicative 
constant k, the thermal conductivity coefficient) equal to the work 
done per cycle in oscillating the surface. The latter result can be 
independently verified by computing the work done per cycle from the 
shear stress on the plate. 


RESULTS AND DISCUSSION 


Now that the temperature distribution and its gradient at the sur- 
face are known, the effect of the surface oscillations on the heat 
transfer can be obtained by comparing any appropriate parameters with 
the corresponding ones for the case of pure conduction to or from a 
stationary surface subject to the same initial temperature difference 
TE. = T) The specific formulation of the parameters depends on the 


particular configuration considered. With this in mind, the parameters 
considered herein will fall into two categories: total parameters, that 
is, those pertaining to configurations in which the oscillations are 
inherent; and net parameters for configurations in which the work done 
in oscillating the surface is assessed to the system. For example, the 
average total heat-transfer coefficient as given essentially by equation 
(14) can be compared with that for conduction to or from a stationary 
surface at an initial temperature difference (Ta ~ T): The temperature 
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distribution for the latter case is given by the sum of the first and 
last terms of equation (13) (see ref. 5, p. 34) and the associated aver- 
age temperature gradient is given by merely the last term of equation 
(14). Therefore, the ratio of these average total heat-transfer coef- 
ficients R, is given by 


U“Pr z Ux ( 16 ) 


1 =l- 


It can thus be seen that the oscillations could lead to higher average 
total heat-transfer coefficients depending on the sign of (Te - T) and 
on the relative magnitude of the two terms in the bracketed part of equa- 
tion (l6); for example, in the limiting case just after the temperature 


Ty is imposed on the surface, that is, tọ + 0, equation (16) reduces to 


Ra #1 a j a) (17) 


10 ~ a 
40 (Tan T) 


so that for the case where the ambient temperature is greater than the 
surface temperature (To >t), and Pr<x", then Rj 71; that is, the 
oscillations increase the average total heat-transfer coefficient. Con- 
versely if Ty? Tes Rio would be larger than unity if Pr- xÊ. 


It should be reiterated that the parameters for comparison depend 
upon the configuration considered. Therefore, if in the specific appli- 
cation the net energy obtained from this system (rather than the total 
as before) is of interest then the total heat transfer mist be properly 
assessed to take account of the work done in oscillating the surface. 
Hence, the net heat transfer for this case (as represented by the dif- 
ference between eqs. (14) and (15)) will be compared with that for con- 
duction to or from a stationary surface. The ratio Rp of the average 
net heat transfer for the oscillating surface to that for a stationary 
surface can therefore be written as 


or (55 
ho VR yy e 
27 (=) 40, (Ta T) 
w sta 


(18) 
oy 


It can be seen from equation (18) that the average net heat-transfer 
coefficient for the oscillating surface would be larger than that for 
the stationary surface (that is, Ro >1) only when Ty? Te and it would 
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be lower when I,<T,. For a given temperature difference, therefore, 
the ratio of the average net heat-transfer coefficients, as given by 
equation (18), can be altered by changing the amplitude of the 
oscillations. 


If the primary interest is not in the net heat transferred across 
the surface (ee it would be, for example, for regenerators) but is in 
the thermal state of the fluid itself, the effect of the oscillations 
can be studied by comparing the average static temperature defect de- 
fined as 


tot £ = 


AT = a | ni (Ta - Tay at (19) 


for the oscillating surface to that for the stationary surface. The 
ratio of these defects is given by 


3 sy" 


“Tose U“Pr SnU“Pr Prg 
a cae 4c,(Ta-Ty)  8ne,(Z,-1,) | 2 (to + zf - (ty i 


For the limiting case, tọ +0, equation (19) reduces to (20) 


2 
Rgg m2 - e ha- 5 fe] (21) 
pe "Ww 


It can thus be seen that unless Pr>64/9 the oscillations will not be 
advantageous for lowering the fluid temperature for the case when 


ie ee 


A more meaningful parameter to examine would perhaps be the ratio 
of the average total-temperature defect with the oscillations to that 
for conduction to the stationary surface. The average total-temperature 
defect is defined as in equation (19) except that T is replaced by the 
total temperature Tų, which is given by 


u 
U 
p 
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in this way the thermal condition of the fluid independent of its motion 
is indicated. This ratio is given by 


oe U“Pr 14 G1 - Pr) [Tx 1 
tc (To 2 Ta) an anPr 5 5 


2 2 


(+ zN- (te) 
(22) 


Thus, the ratio R, could be greater than unity dependent on the rele- 


tive orders of magnitude of the terms in the braces in equation (22) 
and in the limiting case as to +0. 


U“Pr 3 (1 - Pr 
MoT i- Goa R) TE (23) 


Thus, in this case if T,>T, the fluid would be at a lower relative 
temperature with the oscillations (Ryg>1) only if Pr>9. 


An analogous discussion, of course, follows for increasing the 
fluid temperature. 


CONCLUDING REMARKS 


The results of this preliminary analysis of the tempereture distri- 
bution over conducting oscillating surfaces have shown that the heat 
transfer associated with the oscillating surface can be significantly 
different from that for conduction between the same initial temperature 
difference and a stationary surface. In addition to the heat transfer, 
the thermal state of the fluid was studied by means of average static- 
or total-temperature defects and it was found that the surface oscil- 
lations can alter the thermal state of the fluid. In each case, the 
range of parametric values for which the surface oscillations are bene- 
ficial were indicated. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 19, 1954 
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APPENDIX - SYMBOLS 


The following symbols are used in this report: 


Sy 


erf (8) 


k 


Pr 


specific heat at constant pressure 


B 
-5 [ exp(-r)ar 


thermal conductivity coefficient 


frequency of oscillations 


Prandtl number 


ratio of 


ratio of 
to +0 


ratio of 
ratio of 


ratio of 
to + 0 


ratio of 


ratio of 


average 


average 


average 
average 


average 


average 


average 


total heat-transfer coefficients 


total heat-transfer coefficients evaluated as 


net heat-transfer coefficients 
static=-temperature defects 


static-temperature defects evaluated as 


total~temperature defects 


total-temperature defects evaluated as 


static temperature 


equilibrium (or steady-state) static temperature 


total temperature 


surface temperature 


surface temperature for steady state 


static temperature in ambient fluid 


temperature defect 
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t time 

to time at the start of an arbitrary cycle 

U amplitude of oscillations 

U velocity parallel to surface 

y coordinate normal to surface 

œ thermal diffusivity 

T) von Mises coordinate normal to surface 

u absolute viscosity 

v kinematic viscosity 

p density 

Subscripts: 

c complementary solution 

osc oscillating surface 

P particular solution 

sta stationary surface 
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